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The determination of scalar leptons and gauginos masses is an important part of the 
program of spectroscopic studies of Supersymmetry at a high energy linear collider. In 
this talk we present results of a study of the processes: e+e~ — > e+; — > e+e~ Xi X?, 
e+e' p,+ Ji R e+e" Xi Xi, e +e~ &t e~ -»• e+ e~ xl X% and e+e ~ v e v e 
e+ e~ Xi XT m a Supersymmetric scenario at 3 TeV at CLIC. We present the per- 
formances on the lepton energy resolution and report the expected accuracies on the 
production cross sections and on the Sr, p,R, u e , xt an d X? mass determination. 

1 Introduction 

One of the main objectives of linear collider experiments is the precision spectroscopy of new 
particles predicted in theories of physics beyond the Standard Model (SM), such as Super- 
symmetry (SUSY). In this talk, we discuss the production of the supersymmetric partners of 
the muon, electron and neutrino within the so-called constrained Minimal Supersymmetric 
extension of the SM (cMSSM). The parameters chosen are such that the lightest neutralino 
has a mass of 340 GcV, the charginos and heavier ncutralinos have masses in the range 
643 to 917 GeV, the right handed selectron and smuon have a mass of 1010.8 GcV, the left 
handed selectron and the sneutrino have masses of 1110.4 GeV and 1097.2 GeV respectively. 
Smuons are produced in pair through s-chanrrel 7/Z exchange, selectrons and sneutrinos are 
pair produced through s-channel 7/Z exchange or i-channel Xi an d xt exchange respec- 
tively. The branching ratio £^ — > x? is ~ 100 %, and the branching ratios: 6l —J- e~ Xii 
6l — > e~ x!>j an d v c — > e - xt are 16%, 29% and 56% respectively. The cross sections, the 
decay channels and the cross sections times the branching ratio of the signal processes under 
study are given in Table [T] 
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Table 1: Signal processes, cross sections (cr), decay modes, cross sections times branching 
ratio (cr x Br) and cross sections times branching ratio into two electrons and four quarks, 
(cr x Br(ee4Q) ). 
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For the processes e+e" — > t^l^ each t+ decays into a SM lepton and a Xi'i the exper- 
imental signature is two oppositely charged leptons plus missing energy. For the processes 
e+ e" — > e£ e^ — > e+ e~ %2 X2 ~~ ^ e+ e ~ Xi Xi ' l ° ^° an( i e+e~ — > ;> e v e — > e + e _ x+ xT ~~ ^ 
e+ e~ X? Xi W 7- ! the signature is a pair e+e - , four jets and missing energy. The mea- 
surement of the lepton energy distributions of these four processes allows to determine their 
production cross sections and the en, jjR, v e , xt and Xi masses. The aim of this study 
is study, done for the CLIC CDR [T] is to assess the accuracy, which could be obtained on 
the mass measurements and on the production cross sections and characterize the detector 
performances, namely lepton energy resolution. 

2 Event Simulation and Reconstruction 

SUSY signal events and SM background events arc generated using Whizard 1.94 [5] which 
was interfaced to Pythia 6.4 [3] for fragmentation and hadronization. The simulation 
is performed using the GEANT4-based [4] Mokka program [5] with the CLIC_ILD_CDR 
detector geometry [I] , which is based on the ILD detector concept [6] being developed for the 
ILC. The physics backgrounds simulated for this study are listed in Table [2] The centre-of- 



Process 


Decay mode 


a x Br 


a x Br 






fb 


fb 






no cuts 


cuts 


e+e - — ^ 


H+lT 


81.9 


0.6 


e+e - — > \i+ v e \jT v e 


(j+fj,- 


65.6 


3.5 


e+e~ [i+v^v^ 


fi+fj,- 


6.2 


2.2 


e+e- -> W+uW'v 




92.6 


2.4 


e+e" -> Z°vZ°v 




40.5 


0.002 


e+e- All SUSY - (Ar.Mr) 


^+^- 


0.31 


0.31 


e+e- — > e+e~ 


e+e- 


6226.1 


77.1 


e+e- — > e+v e e~v & 


e+e~ 


179.3 


91.1 


e+e- -> W+uW~v 


e+e~ 


92.6 


2.4 


e+e~ -> Z^Z ^ 


e+e- 


40.5 


0.002 


e+e- -> All SUSY - (e+e^) 


e+e- 


1.04 


1.04 


e+e" -> W+W"Z U 


e+erW+W- 


1.35 


0.61 


e+e" -> Z°Z°Z° 


e+e-Z°Z° 


0.045 


0.023 


e+e - — > SUSY — (e + eL and z> c z> c ) 


e+e" (WW, h°h°, Z°Z°) 


0.77 


0.12 



Table 2: Background processes, decay modes and cross sections times branching ration, 
a x Br, without and with preselection cuts. 



mass energy spread coming from the momentum spread in the linac and the beamstrahlung is 
included using the GuineaPig [7] beam simulation for the CDR accelerator parameters [Hi- 
lt is used as input of WHIZARD in which initial state radiation and final state radiation 
(FSR) are enabled. An integrated lumnosity of 2000 fb -1 was assumed, corresponding to 
~3.5 years (1 year = 10 7 s) of run at the nominal CLIC luminosity of 5.9xl0 34 cm" 2 s _1 . 
Events are subsequently reconstructed using the Marlin reconstruction program [9]. The 
track momenta and calorimeter data are input to the PandoraPFA algorithm [10] which 
performs particle identification and returns the best estimate for the momentum and energy 
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of the particles. The muon and electron average identification efficiency are 99% and 96% 
respectively [T|. After reconstruction preselection cuts are applied on the signal samples. 
Table [3] shows the reconstruction efficiencies, £r, for the signal processes. 
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Table 3: reconstruction efficiency, cr without and with 77 overlay and selection efficiencies 
es for the different signal processes. The statistical error on these efficiencies is ~ 1% 



For the process e + e~ — > AhAh there is an inefficiency of about 2.5%; 2.0% is due to the 
cut on the lepton angle and 0.5% is coming from muon misidentification. For the process 
e + e~ — > ejje^ there is an inefficiency of 5.4%; 4.1% is due to the cut on the lepton angle 
and 1.3% is coming from electron misidentification. For the processes e + e _ — > e£e£ and 
e + e~ —¥ v e v e , the parton topology signature required is two leptons and four quarks. After 
reconstruction of all the particles in the event, the jet finder program Fast jet [11] is used 
to reconstruct jets. The jet algorithm used is the inclusive anti-kt method [12] requiring a 
minimum jet energy of 20GcV. An event is retained if 6 jets are found and if two of the 
jets are identified as isolated leptons. Table [3] shows the reconstruction efficiencies of both 
processes, £r is the number of reconstructed 6 jet events, with two leptons, divided by the 
number of generated events with two leptons and four quarks. The energy of the lepton 
is reconstructed from the momentum of the charged particle track corrected for final state 
radiation and brcmsstrahlung; the energy of photons or e + e~ pairs from conversions within 
a cone of 20° around the reconstructed lepton direction is added to the track momentum. 
The lepton energy resolution is characterized using: SE / E^ lue , where SE = -©True — E-n cco , 
©True is the lepton energy at generator level, before final state radiation or brcmsstrahlung 
and ©Reco is the reconstructed lepton energy with photon radiation corrections. Figure [T] 
shows the lepton energy resolutions, for the four processes, the distributions are fitted using 
two Gaussian functions. For the muon final state process, the energy resolution of the peak 
distribution is 5E/E^ UC = 1.5 x 10~ 5 , the energy resolution of the second distribution 
is 5E/E^ mc = 4.9 x 10 -5 and the number of events in the tails is small, 4.1%. For 
the electron final state processes, the energy resolution of the peak distribution is also 
<5_E/i? t 2 mcTi . uc = 1.5 x 10 -5 , but the energy resolution of the second Gaussian is worse, due 
to significant photon radiation, SE/E^, = 8.1 x 10 -5 and the number of events in the tails 
is ~ 30 %. Despite the presence of four jets, the electron energy resolution for the processes 
with two electrons and four jets, Figure [1] (c) and (d), is similar to the one of the process 
with only two electrons in the final state, Figured] (b). To investigate the effect of beam 
induced backgrounds, the reconstruction software is run, overlaying particles produced by 
77 — > hadrons interactions pQ . A sample of 77 — > hadrons events was generated with Pythia 
and simulated. For each physics event the equivalent of 60 bunch crossings of 77 — > hadrons 
events are selected. The detector hits from these events are merged with those from the 
physics event before the reconstruction. A time window of 10 nsec on the detector integration 
time is applied for all detectors, except for the HCAL barrel for which the window is 100 
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Figure 1: Lepton energy resolution, for the processes: e + e — > fi>%fi R (a), e + e — > e^e 
(b), e + e" — > z> e i> e (c) and e + e" — > eje^ (d) 



nsec. After particle reconstruction timing cuts in the range of 1 to 3 nsec are applied in order 
to reduce the number of particles coming from 77 — > hadrons interactions. The cut values 
vary according to the particle type, (photon, neutral hadron, charged particle), the detector 
region, (central, forward) and the Pt of the particle. Table [3] shows the reconstruction 
efficiencies, without overlay and with overlay, after application of the selection cuts. The 
time selection cuts preserve the energy resolution, but induce an inefficiency of 6% for 
the dielectron final state processes, 0.5% is coming from the current pattern recognition and 
tracking software, the rest is due to the increase of lepton misidentification. For the processes 
e + e" — > /U^/ifl ail d e + e~ — > ^r,^ri n0 timing cuts are applied. The energy resolution of 
these processes is not affected by the 77 — > hadrons interactions because the high cut of 
4 GeV, applied to reduce the SM physics background, removes all the hadrons coming from 
these interactions; it induces a reconstruction inefficiency of 1.0% and 4.6% respectively, see 
Tabled 
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3 Event Selection and Mass Determination 



To tell apart signal events from background events the following set of discriminating vari- 
ables is used: dilcpton energy E(L1) + E(L2), vector sum Pt(L1) + Pt(L2), algebric sum 
Pt(L1) +Pt(L2), dilcpton invariant mass M(L1,L2), dilepton velocity /3(L1,L2), angle of 
the dilepton missing momentum vector cos#(Ll, L2), dilepton acolinearity 7r — 9 2 — 0\, dilep- 
ton acoplanarity tt — </> 2 — 4>\ and energy imbalance A = |E(L1) — E(L2)|/|E(L1) + E(L2)| 
where LI and L2 are the two leptons. Histograms of the discriminating variables are built 
for signal and background events. The events are weighted such that the data samples cor- 
respond to the same integrated luminosity. From these histograms, signal and background 
event probability density functions are computed and combined into a total probability 
classifier using the multivariate analysis toolkit, TMVA [13]. The signal and background 
samples are split into two equal event samples called "Monte Carlo" and "Data". The 
Monte Carlo sample is used to train the classifier which ranks events to be signal or back- 
ground like. The method is then applied to the so called Data sample, for each event a 
total probability is computed and a cut is applied to tell apart signal from background. The 
cut value is chosen to optimise the significance Ns/a/(Ns + N B ) versus the signal efficiency; 
Ng and Nb are the number of signal and background events. The selection efficiencies, e$ 
are shown in Table [3] After selection and background subtraction, the slepton, neutralino 
or chargino masses are extracted from the position of the kinematic edges of the lepton 
energy distribution, a technique first proposed for squarks [14], then extensively applied to 
sleptons [To] . 

{E H -E L fV /2 / 2(E H +E L )^ 1/2 



= (Eh + El) 2 J ^ ™*? ° r ^ = m i- ^ ^ ) « 

The masses depend on the centre-of mass energy \fs/2 and on the kinematic edges values 
El^h, therefore the accuracy on the masses relies on the measurement of the shape of the 
luminosity spectrum and on the lepton energy resolution. The masses arc determined using 
a 2-parameter, m^ ± and m^a, x 2 fit to the reconstructed energy distribution. The fit is 
performed with the Minuit minimization package }16| . The lepton energy spectrum is a 
uniform distribution with end points fixed by the slepton and neutralino masses. For each 
event a random value of y/s is generated taking into account the beamstrahlung and ISR 
effects; the lepton energy resolution is included using the energy resolution functions shown 
in Figure [1] The process cross section is obtained from the integral of the momentum 
distribution. Table 0] shows the values of the measured sleptons cross sections, sleptons 
masses and gauginos masses. Figure [2] (a) and (b) show the lepton energy distributions 
and fit results for the processes, e + e~ — > MjjAfl an d e + e~ — > v e v e respectively. For the 
process e+e~ — s- — s- e+ e~ x° X2: the cross section is determined from the fit to Di-jct 
invariant mass shown in Figure [2] (c). 
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Table 4: Cross sections values, sleptons and gauginos masses and statistical accuracies. 
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Figure 2: Lcpton energy spectrum and fit results, for the processes: (a) e+e - —> fl^p.^, (b) 
e+e - — > v e v e: (c) Di-jet invariant for e+e - — > v e v e and e+e - — > eje£ 
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